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Abstract. Particle number fluctuations and correlations in nucleus-nucleus collisions 
at SPS and RHIC energies have been studied within the Hadron-String-Dynamics 
(HSD) transport approach. Event-by-event fluctuations of pion-to-kaon, proton-to- 
pion and kaon-to-proton number ratios are calculated for the samples of most central 
collision events and compared with the available experimental data. It has been found 
that the HSD model can qualitatively reproduce the measured excitation function for 
the K/ir ratio fluctuations in central Au+Au (or Pb+Pb) collisions from low SPS up 
to top RHIC energies. These predictions impose a challenge for future experiments. 
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1. Introduction 

The measurement of the fluctuations in the kaon to pion ratio by the NA49 Collaboration 
was the first event-by-event analysis in nucleus- nucleus collisions [1]. It was suggested 
that this ratio might allow to distinguish events with enhanced strangeness production 
attributed to the QGP phase [2]. Nowadays, the excitation function for this observable 
is available in a wide range of energies: from the NA49 collaboration [3] in Pb+Pb 
collisions at the CERN SPS and from the STAR collaboration [4,5] in Au+Au collisions 
at RHIC. Results from NA49 show an enhancement of fluctuations in the kaon to pion 
multiplicity ratio for low energies which may be a signal of a deconfinement phase 
transition. On the other hand there is no such an enhancement for the proton to pion 
ratio fluctuations. 

This proceeding presents the results of a systematic study of K/tt, K/p and p/w 
ratio fluctuations (K = K + + K~, n = 7r + + 7r~, and p means p + p) based on 
the HSD transport model [6]. For more details we refer readers to [7,8]. Event-by- 
event fluctuations for charged hadron multiplicities also have been studied in the HSD 
transport approach [9]. 

Let's introduce some notations. The deviation ANa from the average number (Na) 
of the particle species A is defined by N A = (Na) + ANa, while the covariance for species 
A and B is: 

A (Na, N b ) = (ANaAN b ) = (N A N B ) - (N A ) (N B ) , (1) 

the scaled variance and the correlation coefficient 

A(N A ,N A ) (N 2 a )-(N a ) 2 (ANaANb) 



UA = (N A ) ~ (N A ) ' PAB 



{(AN A ) 2 ) ((AN B ) 2 ) 
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• (2) 



The fluctuations of the ratio Rab = Na/Nb will be characterized by a 2 [2,10], which 
to the second order in ANa/ (Na) and ANb/ (Nb) can be written as: 



2 _ ((ARab) ) „ uj A ujb „ 

° = (Rab) 2 = W) Wb~) " 2PAB 

2. Mixed Events Procedure 



ua^b 
(N a )(N b ) 
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(3) 



The experimental data for Na/Nb fluctuations are usually presented in terms of the so 
called dynamical fluctuations [11] 

a dyn = signed - a mix ) a - a mix , (4) 

where a 2 is defined by (j3J), and a^ix corresponds to the following mixed events procedure. 
One takes a large number of nucleus-nucleus collision events, and measures the numbers 
of Na and Nb in each event. Then all A and B particles from all events are combined 
into one set. A construction of mixed events is done like the following: One fixes a 
random number iV = A^^ + iV^ according to the experimental probability distribution 
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P(N), takes randomly N particles (A and/or B) from the whole set, fixes the values of 
Na and Ng, and returns these N particles into the set. This is the mixed event number 
one. Then one constructs event number 2, number 3, etc. 

Note that the number of events is much larger than the number of hadrons, N, 
in any single event. Therefore, the probabilities pa and ps = 1 — Pa to take the A 
and B species from the whole set can be considered as constant values during the event 
construction. Another consequence of a large number of events is the fact that all A and 
B particles in any constructed mixed event most probably belong to different physical 
events of nucleus-nucleus collisions. Therefore, the correlations between Nb and Na 
numbers in a physical event are expected to be destroyed in a mixed event. This is the 
main purpose of the mixed events construction. 

Calculating the Na/Nb fluctuations for mixed events one gets like in [7]: 

2 _ 1 1 

The mixed event procedure gives the same a 2 for Na/Nb fluctuations as in the GCE 
formulation for an ideal Boltzmann gas, i.e. u>a = ojb = 1 and pab = in d3J). Note 
that non-trivial fluctuations of Na and Nb (w™J ^ 1) as well as non-zero correlations 
pmix ma y ex j S £ j n the mixed events procedure, but they are canceled out exactly in cx^.. 



3. HSD Results and Comparison to the Data 

The HSD results as well as statistical model calculations for ratio fluctuations in full 
acceptance are presented in [7, 8] A comparison of the SM results for fluctuations in 
different ensembles with the data looks problematic at present due to the difficulties 
in implementing the experimental acceptance and centrality selection which, however, 
can be taken into account in the transport approach. In order to compare the HSD 
calculations with the measured data, the experimental cuts are applied for the simulated 
set of HSD events. In Fig. [I] the HSD results of a^ yn for the K/ir, p/n and K/p ratios 
are shown in comparison with the experimental data by the NA49 Collaboration at the 
SPS [3] and the preliminary data of the STAR Collaboration at RHIC [4,5,12,13]. The 
available results of UrQMD calculations (from Refs. [3, 14, 15]) are also shown by the 
dashed lines. 

For the SPS energies we use the NA49 acceptance tables from Ref. [3]. For the 
RHIC energies we use the following cuts: in pseudorapidity, \rj\ < 1, and in transverse 
momentum, 0.2 < pr < 0.6 GeV/c for kaons and pions and 0.4 < pt < 1 GeV/c for 
protons [4,5,12,13]. Note, that HSD results presented in Fig. [1] correspond to the 
centrality selection as in the experiment: the NA49 data correspond to the 3.5% most 
central collisions selected via veto calorimeter, whereas in the STAR experiment the 5% 
most central events with the highest multiplicities in the pseudorapidity range \rj\ < 0.5 
have been selected. 

One sees that the UrQMD model gives practically a constant cr^j, which is by 
about 40% smaller than the results from HSD at the lowest SPS energy. This difference 
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Figure 1. The HSD results for the excitation function in adyn for the K/n, K/p, 
p/t: within the experimental acceptance (solid line) in comparison to the experimental 
data measured by the NA49 Collaboration at SPS [3] and by the STAR Collaboration 
at RHIC [4,5,12,13]. The UrQMD calculations are shown by dotted lines. 



between the two transport models may be probably attributed to different realizations 
of the string and resonance dynamics in HSD and UrQMD: in UrQMD the strings 
decay first to heavy baryonic and mesonic resonances which only later on decay to 'light' 
hadrons such as kaons and pions. In HSD the strings dominantly decay directly to 'light' 
hadrons (from the pseudoscalar meson octet) or the vector mesons p, uo and K* (or the 
baryon octet and decouplet in case of baryon number ±1). Such a 'non-equilibrated' 
string dynamics may lead to stronger fluctuations of the Kj tt ratio. 

At the SPS energies the HSD simulations lead to negative values of Odyn for the 
proton to pion ratio. This is in agreement with the NA49 data in Pb+Pb collisions. 
On the other hand HSD gives large positive values of o v ^ yn at RHIC energies which 
strongly overestimate the preliminary STAR data for Au+Au collisions [12]. For a^ y p n 
only preliminary STAR data in Au+Au collisions are available [13] which demonstrate 
a qualitative agreement with the HSD results (Fig. CQ. The HSD results for a^ yn show 
a weak energy dependence in both SPS and RHIC energy regions. A peculiar feature is, 
however, a strong 'jump' between the SPS and RHIC values, seen in the middle panel 
of Fig. [TJ in the HSD calculations which is caused by the different acceptances in the 
SPS and RHIC measurements. 

The influence of the experimental acceptance is clearly seen at 160 A GeV where a 
switch from the NA49 to the STAR acceptance leads to the jump in a d y n by 3% - middle 
panel of Fig. [TJ On the other hand, our calculations for Pb+Pb (3.5% central) and for 
Au+Au (5% central) collisions - performed within the NA49 acceptance for both cases 
at 160 A GeV - shows a very week sensitivity of a^ yn on the actual choice of the collision 
system and centrality - cf. the coincident open circle and triangle at 160 A GeV in the 
middle panel of Fig. [TJ 
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4. Summary 

The event-by-event multiplicity fluctuations of K/n, K/p and p/ir ratio fluctuations in 
central Au+Au (or Pb+Pb) collisions from low SPS up to top RHIC energies have been 
studied within the HSD transport approach. 

It has been found that the HSD model can qualitatively reproduce the measured 
excitation function for the K/tt ratio fluctuations in central Au+Au (or Pb+Pb) 
collisions from low SPS up to top RHIC energies. Accounting for the experimental 
acceptance as well as the centrality selection has a relatively small influence on o~dyn and 
does not change the shape of the ad yn excitation function. 

The HSD results for a v ^ yn appear to be close to the NA49 data at the SPS. The data 
for a^y n in Pb+Pb collisions at the SPS energies will be available soon and allow for 
further insight. A comparison of the HSD results with preliminary STAR data in Au+Au 
collisions at RHIC energies are not fully conclusive: o~d yn from HSD calculations is 
approximately in agreement with data [13] for the kaon to proton ratio, but overestimate 
the experimental results [12] for the proton to pion ratio. New data on event-by-event 
fluctuations in Au+Au at RHIC energies will help to clarify the situation. 
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